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Abstract 

A general formula is derived for the finite volume dependence of vacuum expec- 
tation values analogous to Liischer's formula for the masses. The result involves no 
integrals over kinematic quantities and depends only on the matrix element between 
pions at zero momentum transfer thus presenting a new way to calculate the latter, 
i.e. pion sigma terms. 

The full order correction to the vacuum condensate (qq) is evaluated and 
compared with the result from the Liischer formula. Due to the size of the result 
no conclusion about the accuracy of the Liischer formula can be drawn. 



1 Introduction 



Quantum Chromo Dynamics (QCD) at low energy remains a difficult problem. One of 
the ways to deal with this problem is to numerically evaluate the functional integral of 
QCD. This approach, known as lattice QCD, has now reached the stage where realistic 
calculations with fairly light quark masses are now possible. One side effect of this is 
that finite volume corrections are becoming more important. Luckily in many cases these 
corrections can be evaluated analytically using Chiral Perturbation Theory (ChPT) PIE]- 
The application of ChPT to finite volume was started by Gasser and Leutwyler |3]. A 
review of recent work in this area can be found in jl]. Note that ChPT is applicable to 
finite volume as soon as the typical momenta that are relevant are small enough. This 
imposes a size restriction on the volume as 



Here F^^ is the pion decay constant and L is the linear size of the volume. This paper is 
concerned with the p-regime. This is the regime where the volume is large enough such 
that the zero momentum fluctuations of the meson fields can be treated perturbatively. 
This is the regime with in addition the requirement that 



These finite volume corrections have been evaluated for many quantities up to one-loop 
order. This is the order where the ffist nontrivial dependence on the volume shows up. 
One purpose of this paper is to calculate the full two-loop finite volume corrections to the 
vacuum condensate {qq). This is the one of the calculations of finite volume effects to this 
order^ The vacuum condensate at finite volume has been studied at one- loop in Ref. [B]. 

An alternative approach to finite volume corrections was introduced by Liischer where 
the leading part of the finite volume corrections was derived to all orders in perturbation 
theory for the mass in terms of a scattering amplitude [71 |H]. This was extended to the 
finite volume corrections for the pion decay constant in [Oj. The other purpose of this paper 
is to extend the Liischer formula also to vacuum expectation values. This will in general 
connect the finite volume corrections of an operator to the zero-momentum transfer matrix 
element of that operator between pion states as shown in Eq. (0). This allows for new ways 
to calculate sigma terms from the finite volume variation of vacuum condensates. 

Note that all our formulas are for the case of an infinite extension in the time direction 
but a finite volume in the three spatial directions. The formulas can be easily extended 
to a small fourth direction as well by replacing the integrals with the expressions for that 
case. 

^Ref. [ni with another calculation appeared essentially simultaneously. There the pion mass at finite 
volume was studied at two-loop order. 



F^L>1. 



(1) 



mlF^V » 1 . 



(2) 
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2 A Liischer Formula for the Vacuum Condensate 



It turns out to be straightforward to extend Liischer's formula for the mass to the case of 
the vacuum condensate. In fact the formula has an even simpler structure than the one 
for the mass or the decay constant |9^. The underlying observation of Liischer's method 
is that the leading finite volume corrections come from one of the propagators feeling the 
finite volume effects only. We write the finite volume propagator Gv{q) as 

Gy(g) = Ee-^'''""^^G'oo(g). (3) 

n 

Here the sum is over a vector of integers and Lj is the length of the volume in the i-th 
direction. The term with n = o gives the infinite volume case. Liischer then kept only the 
the term with one of the = ±1 since this is parametrically leading. As suggested in pH] 
one can also easily keep the entire series, thus keeping part of the subleading corrections. 
The remainder of the Liischer's analysis goes by taking the component of q parallel to ra, 
Qn and distorting the integration contour along that direction to q-a ~^ Qn — is. In the limit 
of s oo the contribution from that integral vanishes and only the parts coming from 
the singularities encountered while deforming the contour remain. These happen when 
relevant propagators can go on-shell. A pedagogical introduction can be found in 8j. 

The difficulty in proving this is that it needs to be done to all orders in perturbation 
theory, this was done in 0. Here we only sketch the lines of reasoning. Of the three types 
of contributions shown in Fig. 4 of [7j only one is relevant here and is shown in Fig. ^b). 
The case of operators with contributions of the type shown in diagram (a) is not relevant 
for ChPT, one can always add the parity-conjugate operator to remove the contribution 
from diagram (a). 

One difference with the case of the mass or decay constant here is that when the 
singularity is encountered, there is no freedom left in the relevant matrix element and the 
integral over the momentum of the propagator can be fully done. In the case of the mass 
and decay constant there is an external momentum available, Pexti this leaves after putting 
the deformed q on-shell a freedom in q ■ pext which results in the final integration over v. 
Here there is no external momentum and hence there is no such freedom left. 

The formula after putting everything on-shell becomes 

{0)v - {OU = -11-^2 r /"^f^' e-V»-K+^-)^-(0|O|^) , (4) 

for a real boson with infinite volume mass uiq and the matrix element on the right hand 
side should be taken at zero momentum transfer. The integral can now be done explicitly 
in terms of the generalized Bessel function Ki. It also only depends on k = which is 
also integer. The number of times in the sum that = k we call x{k). We obtain 

{o)v-{ou = - E ^2^M((k)Mo\<P) , (5) 

fc=l,oo V^V ) 
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(a) 



(b) 



Figure 1: The two types of diagrams giving the leading finite volume corrections 



with ((k) = ^/km^L. Note that the above formula is for the case of one real scalar. The 
multiplicity factors for complex scalars can be trivially taken into account. Note that we 
have left the sum over all modes in as suggested for the mass in pjGj . 

In particular for the case of {qq) the relevant matrix element is the sigma term. The 
finite volume corrections to the vacuum condensate thus are another option to calculate the 
pion sigma term. This can then be compared with the direct calculation of the sigma term 
via the matrix element (vr|gg|7r) or via the Feynman-Hellman theorem from dm^/dniq. 



3 The finite volume vacuum condensate at two-loops 

The vacuum condensate at two- loop was calculated in ^2] , here we repeat that calculation 
taking into account the finite volume effects. The calculation in terms of the lowest order 
meson masses is straightforward and proceeds exactly as in J2] • The details of calculation 
of two loops in ChPT can be found in The diagrams that contribute are shown in 

Fig.ia 

In infinite volume the loop diagrams contain the integrals 

1 f d'^p 1 



B{m') = - f (6) 
These are expanded in terms of e = (4 — (i) /2 as 

= ^f--7ir + log(47r) + l))+]4K) + eA^(m2) + 0(e2), 
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Figure 2: The diagrams up to order for (gg). The hues are meson propagators and the 
vertices are: o a. insertion of gg, ® a. p^ insertion of gg, a insertion of gg, • a 
vertex and x a vertex. 

= -i-^Q-^^ + log(47r) + l))+5K) + e5^(m2) + 0(e2), (7) 

with a similar expansion for B{m?). At finite volume the integrals over momenta get 
replaced by sums over the finite possible momenta. In this paper we only keep three of the 
four dimensions at the same finite length L. Their principle of evaluation can be found 
in Refs. jHl and explicit expressions are at finite volume. After renormalization the 
divergent parts cancel and the finite parts become subtraction-scale dependent. The 
needed explicit expressions are 

A{m^) = -i^logf!^!-^ y x{k)^K,{Xk) 

Bim^) - -T^l°^(f)-T^ + T^ E^-MA.) (8) 

with Afc = ^/k mL . The functions Ki and K2 are modified Bessel functions and the integer 
quantity x{k) is the number of times the sum of squares k"^ = kf + k^ + k1 is equal to k 
when ^1,^2 5 ^3 are varied over all positive and negative integers. For large L the Bessel 
functions lead to an exponential fall-off with the finite size. 

We have checked explicitly that at two-loop order the contributions containing A"^ and 
cancel, so we do not need to evaluate these at finite volume. That this cancellation 
happens follows from the cancellations of nonlocal divergences but in the calculation of ^2] 
this calculation was not explicitly checked. We have also reperformed the full calculation 
since in ^21 the relation Aim?) = m? {B{rn?') + l/(167r^)j was used, which is no longer 
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true at finite volume. In Ref. the result was expressed in terms of the physical mass 
and decay constant. This can be done here as well but since we must then distinguish 
between the integrals at finite volume and those at infinite volume the formula is shorter 
if we leave meson masses and decay constants in their unrenormalized form. 

The result, in three flavour ChPT, uses the low-energy constants (LECs) L[ from 
Ref. [2j and the from [14j. The latter drop out when comparing finite and infinite 
volume corrections and the values of the L[ only contribute to that difference at NNLO. 
We also suppress the dependence of all quantities. The lowest order masses for the pions, 
kaons and eta we denote by Xn, Xk and x?? respectively and in term of the average up and 
down quark-mass m and strange quark-mass rris they are 

Xn = 2Bom, Xk = Bo{m + nis) Xv = '^^oim + 2ms)/3 . (9) 

We define the quantities 

(11 u \ / s s \ 

The calculation gives 

< = A{xn) {s/2) + A{xk) + A{x,) {l /6) + 16 Xn LI + 8 XnLl + 4 XnH; + 32 XK LI, 
vl = A{xk){2)+A{x,){'2/s)+16x.LI + 4{2xk-X.){2LI + H;)+32xkLI, 

< = A(x.)'(-3/8)+A(x.)A(x,)(l/4)+A(x.)5(x.)(-3/4x.) 

+A{xn) Bixv) (l/12 Xn) + A{xk) Mx,) ( - l/s) + Mxk) B{Xn) ( " 2/9 Xk) 
+A{x,y (l/72) + A{x,) B{Xn) (l/4 Xn) + A{x,) B{Xk) ( - 1/3 Xk) 
+A{x,) BiXv) ( - 7/108 X. + 4/27 xk) 

+A{xn) {-36X.LI- 24 x. L; + 72 Xn L^ + 48 x. - 24 Xk L^ + 48 xk L',) 
+A{xk) [~8x.Ll + lQx. LI - 48 xk LI - 16 xk LI + 96 xk + 32 xk L;) 
+A{x,) (4/3 X. LI + 8/9 Xn L; - 8/3 Xn L^ + 64/3 Xn L'j + 16/3 x. L; 

-40/3 Xk LI - 32/9 xk L; + 80/3 xk L; - 64/3 xk L't) 
+B{x.) ( - 24 X. XK LI + 48 Xk LI - 12 xl LI - 12 xl LI + 24 xl + 24 xl Ll) 
+B{xk) {-8x.XkLI + 1Q Xn Xk L; -16x1 LI -8 x\ ^5 + 32 x\ LI + 16 xl ^s) 
+B{xn) ( - 8/9 Xn XK LI + 32/27 Xn Xk L; + 16/9 Xn Xk L^ - 128/9 Xn Xk L\ 

-64/9 Xn XK LI + 4/9 xl LI - 4/27 xl ^5 " 8/9 xl LI + 64/9 xl ^ + 8/3 xl LI 
-32/9 xl LI - 64/27 xl ^5 + 64/9 xl LI + 64/9 xl L'j + 64/9 xl Ll) 
+ 192 Xn Xk C;, + 8XnXK C^, + 48 xl C\., + 80 xl + 48 xl ^^21 
-4 xl CI, + 64 xl CI, + 192 xl C21 , 
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v', = +Aix.) Bixr,) (1/3 X.) + Mxk) Aixv) ( - 2/3) + A(xi^) ( - 8/9 xk) 

+A{x,y (2/9) + ( - 2/3 Xi^) + ( - 7/27 + 16/27 xk 

+A{x.){-24x.Ll + 48xnLl) 

+A{xk) {-IQX.LI + 32 X. LI - 64 xk LI - 32 xk L; + 128 xk L^ + 64 Xi^ i^^) 
+:4(x,) ( - 8/3 X. + 32/9 x. ^5 + 16/3 x. - 128/3 x. - 64/3 x. 
-64/3 XK LI - 128/9 xk LI + 128/3 Xi^ LI + 128/3 Xi^ + 128/3 xk L^) 
+B{xk) ( - 16 X. XK LI + 32 X. Xi^ - 32 xl LI - 16 x^ L; + 64 x^ 
+32x^^8) 

+B{xr,) ( - 32/9 Xn Xk L^ + 128/27 x. Xi^ ^5 + 64/9 Xn Xk L; - 512/9 x. Xk L'7 
-256/9 X. XK LI + 16/9 x^ LI - 16/27 x^ - 32/9 xl L^ + 256/9 xl L'7 
+32/3 x^ L; - 128/9 x^ - 256/27 xl ^5 + 256/9 xl LI + 256/9 xl ^7 
+256/9 x^i^S) 

-192 X. Xk CI, - 64 x. Xk + 192 X. Xk C;, + 48 x^ C^, + 16 x^ C^^o 

+48 xl CI, + 4 xl CI, + 192 x^ CI, + 192 x^ C^o + 192 x\ C2, . (11) 

These results agree analytically with those of Ref . ^2] • In Ref . numerical results were 
presented. Using the formulas above one obtains much smaller numerical corrections at 
NNLO then were obtained there. This effect is mainly due to the rewriting of the I/Fq 
into 1/-F^ and to a lesser extent of rewriting the masses in terms of the physical masses. 
Numerical results are presented in terms of the ratio 



{m )v - { qQ)c 



where we calculate both numerator and denominator to NLO or NNLO in ChPT. As 
input parameters for Fq and the L[ we use the values obtained in fit 10 of Ref. ^H]- In 
addition, we have set H2 = 0. The results for both and Rg are shown in Fig. |21 for 
Xk = (450 MeV)^ and three values of the lowest order pion mass Xn = (100 MeV)^, 
(250 MeV)^ and (450 MeV)^. The finite volume corrections to the strange quark vacuum 
expectation value are always small. The light quark vacuum expectation value can have 
sizable finite volume corrections for the smaller pion mass. 



4 Comparison and conclusions 

The sigma terms in infinite volume are known to two-loop order in ChPT. Either directly 
[TB] or via the derivative of the meson mass to the quark mass ^21 ^\ ■ The sigma terms 
are also known for two- flavour case p!Hl IT^ 120] ■ 
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(a) (b) 

Figure 3: (a) The ratio Ru for three values of the input lowest order pion mass and 
Xk = (450 MeV)2 (b) The same for R,, with xk = (450 MeV)^ and x^ = (100 MeV)^, 
and (450 MeV)^. 



A major motivation for this work was to test the accuracy of the Liischer type of finite 
volume formulas. So, how well do the two approaches compare. 

At one-loop the comparison is rather trivial since at most one-propagator can show 
up in the relevant one-loop diagrams. We can also apply the formula (0) for the different 
species of meson separately and thus construct the exact one loop ChPT formula. At one 
loop the only test possible is thus how fast the sum over k converges and how quickly 
it converges to an exponential. This was already studied in ^T] and the convergence to 
the leading exponential is rather slow while the sum over k converges faster. We find the 
same conclusions. The extended Liischer formula agrees analytically with the full one-loop 
expressions and is thus fully accurate. 

At two loop level the two formulas have a different behaviour. There are diagrams 
now allowing two propagators simultaneously to feel the effect of the finite volume. A 
small complication that needs to be taken into account here is that the Liischer formula 
is with the infinite volume mass at one loop. Thus when changing from the lowest order 
mass to the physical mass in the one loop formulas this needs to be done with the infinite 
volume expressions. But, even after doing this, the corrections are very small. The two- 
loop calculation as plotted in Fig. El is obviously very small. In addition, it is dependent 
on precisely how one defines the one- loop order. E.g., there are ambiguities in using the 
the physical pion decay constant or the lowest order one, and how much one uses the 
Gell-Mann-Okubo relation in the one-loop expression. With these changes the two-loop 
calculation can be changed significantly but remains mostly small. The Liischer formula 
at this order has also small corrections since the sigma terms have small corrections at 
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one-loop. We have not plotted it since it will essentially be on top of the other curves in 
Fig. El 

Both the extended Liischer formula and the full two-loop calculations thus indicate 
small two-loop order corrections. The actual numerical results of the two-loop expressions 
depends strongly on the inherent ambiguities in defining it. So, both approaches give 
comparable results at this order, but we cannot draw conclusions on the accuracy of the 
extended Liischer formula. 

The extended Liischer formula allows also to include effects from even higher orders by 
using the sigma terms at higher orders. This quantity is known to NNLO but its numerical 
value depends strongly on the input parameters chosen (TH]. There might thus be sizable 
effects at higher orders but there need not be. 

In conclusion, we have derived an extended Liischer formula for the finite volume effects 
on the quark vacuum condensate. We have also calculated these effects to two-loop order in 
ChPT. At one-loop order the extended Liischer formula is exactly equal to the full ChPT 
calculation. At two-loop order, the latter includes extra effects, but both approaches 
indicate very small corrections. The difference is within the inherent uncertainty of the 
full ChPT calculation at that order and we thus cannot conclude if the difference is due to 
the inherent uncertainty in the two-loop order result or to the effects not captured by the 
Liischer formula. 
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